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Abstract 

We discuss the possibility of measuring generalized parton distributions in exclusive 
electroproduction of mesons off the nucleon and estimate the uncertainty from pertur- 
batively induced higher-twist corrections. We find that, while the magnitude of the cross 
section changes significantly taking into account twist-four contributions modeled via 
renormalons, the transverse spin asymmetry is weakly sensitive to them and displays the 
precocious scaling. 
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Abstract. 

We discuss the possibility of measuring generalized parton distributions in exclusive electro- 
production of mesons off the nucleon and estimate the uncertainty from perturbatively induced 
higher-twist corrections. We find that, while the magnitude of the cross section changes signifi- 
cantly taking into account twist-four contributions modeled via renormalons, the transverse spin 
asymmetry is weakly sensitive to them and displays the precocious scahng. 

Generalized parton distributions (GPDs) F(x, (^,A^) encode exhaustive information 
on one-particle correlations in the nucleon and thus carry the lore on its wave function 
and the phase structure of the latter. The quantum-mechanical wave function W allows 
to predict expectation values of all observable for a given system. An identical descrip- 
tion is achieved by means of the density matrix p{xijX2) = W*{xi)W{x2). The latter 
can be used in turn to construct the quantum equivalent of the classical phase-space 
distribution, the primary example being known as the Wigner quasi-probability function 
W{k,r) = f ^e^'^^/^p(r— ^x,r+^x). Contrary to its classical counterpart it is not pos- 
itive definite, — a hallmark of the interference. The marginals of W{k, r) acquire how- 
ever the probability interpretation as coordinate density p{r) — J dkW{k,r) — |'F(r)p, 
or equivalently the Fourier transform of the atomic form factor, and momentum-space 
distribution nik) = J ^W{k, r) = with ^{k) = J ^q-'^I^W{x). The W{k, r) 

is an analogue of a Fourier transformed one-dimensional GPD. The impact parameter- 
dependent parton distributions [1], related to GPDs again by a Fourier transform with 
respect to the momentum transfer Aj^, are transparently identified as relativistic nucleon 
Wigner distributions [2]. Thus, the studies of GPDs will shed the light on the phase- 
space distribution of quarks in the proton. GPDs are cleanly probed in deeply virtual 
Compton scattering involving only one hadron, the nucleon, whose structure is unrav- 
eled through electron scattering [3]. The same GPDs enter the amplitude of exclusive 
electroproduction of mesons in the asymptotic regime of large momentum transfer [4]. 
However due to the presence of an extra hadron in the final state and the specifics of per- 
turbative QCD approach to such processes, one has to pose the question of applicability 
of hard gluon exchange mechanism at moderate photon virtualities. The cross section of 
meson photoproduction with longitudinally polarized 7* is, see Fig. 1 (left). 
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FIGURE 1. Kinematics of the exclusive meson electroproduction off the proton in its rest frame (left) 
and leading order perturbative diagrams in hard scattering approach (right) dressed by fermion bubble 
insertions which generate power corrections in the amplitude. 



where q\ = xb = {2pi ■ qi), and M runs over mesons of different flavors. For 
the charged psedoscalar meson P = the decay constant is = 132MeV, while for 
the neutral vector mesons 

yo = pO^ jj^gy ^ 153MeV and fa = 138MeV. The 

parts of the cross sections read for target polarization-independent [5, 6, 7] 
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and target (transverse) polarization-dependent components 
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respectively. The generalized structure function ^ = {Jif.S'.J^.S'} depends on the 
skewness ^ =X'b/{2—xb), the f-channel momentum transfer and resolution scale 
In leading-twist approximation, it is expressed as a convolution of the meson distribution 
amplitude normalized to /q^ Jm^(m) = 1, the quark or gluon GPDF = {H,E,H,E} 
and, correspondingly, the quark or gluon coefficient function T via [4] 

1 rl 

du J ^dx^M{u) {TMiu,x,t,;^^)FM{x,^,A^) + Tg{u,x,t,)Fg{x,t,,A^)} . 



The function Tm to lowest order approximation is given by the one-gluon exchange 
mechanism diplayed in Fig. 1 (right). The studies of higher-order perturbative correc- 
tions to the hard coefficient function in many physical observables have demonstrated 
that ambiguities generated by the perturbative resummation of fermion vacuum polar- 
ization insertions were of the same order of magnitude as available non-perturbative 
estimates of matrix elements of higher-twist operators. The development and sophistica- 
tion of these ideas has led to some evidence that infrared renormalons might reflect the 
magnitude of higher-twist contributions and even their functional dependence on scal- 
ing variables and can thus be used as a rough estimate of power-suppressed effects. On 
the practical side to compute them in the present circumstances, one replaces the tree 
gluon propagator, in the single bubble-chain approximation, see Fig. 1 (right), by [in the 
Landau gauge] 

where Cms = f the MS and Cms = f — 7b + ln4;r in the MS scheme, and b = 

^-Nc — ^TfNf is the first coefficient of the QCD beta-function and = aj(jU^) = 
4n/ (blnji^/A^), where the last equality hold to one-loop order. The functions Fm 
which enter the above structure functions are combinations of ^-flavor quark GPDs 

Fn = Fu-Fd, Fp = QuFu-QdFd, Fco = QuFu + QdFd, U) 

where the quark charges are Qu = ^ and Qd — — ^. For and only the polarized 

F = {H,E} and, correspondingly, unpolarized GPDs F = {H,E} enter the game. The 
quark coefficient function (with resummed renormalon chains) for the M = it'^ has the 
form 
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with Cf = {N^ — \)/2Nc- The coefficient function for neutral vector mesons M = V^, 
Ty, is obtained from this one by setting Qu,Qd^ 1 since the quark charges are included 
into flavor combinations of GPDs. Finally, for completeness we present the leading- 
order gluon coefficient function contributing to neutral vector meson production. 
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If one absorbs the dependence on the momentum fraction into the argument of the 
coupling, as{ju{l ± |-)^^e^'^), one explicitly sees that the end-point regions produce 
divergences. Infrared renormalons are caused by the end-point singularities [Feynman 



mechanism] in exclusive amplitudes [8], see also [9, 10]. This can be viewed as an 
estimate of the ambiguity in the resummation of higher-order perturbative corrections 
or, taken to the extreme, as a model of higher-twist contributions [11]. Convolution of 
the coefficient function with the distribution amplitude generates renormalon poles. For 
the asymptotic distribution amplitude 0asy(M) = 6uu one gets two poles T = 1 and T = 2, 
corresponding to ambiguities on the level of and power corrections. Since the 
latter receives extra contributions from higher order diagrams as well, we use only T = 1 
pole for the estimates of the form of higher-twist corrections. Taking the imaginary part 
(divided by n) arising from the contour deformation around the renormalon poles as a 
measure of their magnitude, we get 

~ ~ A— e^/^ ~ 

, A2; ^ ^^PV(^ , A2; ^ j^^ ^^^^ ^ ^ ^ ^2) ^ 

where = ±l comes from the ambiguity to go around the renormalon pole in the Borel 
plane. Here we have used the one-loop expression for the QCD coupling constant and 

A.(-.«) = 48^{(T3f^-(TTf^)' 

Since the GPD and its first derivative are continuous functions atx = ±^ [12], the second 
integral is well-defined. In the first term of (9) one uses the principal value prescription 
to go around the poles in the Borel plane. For the E we use the pion-pole dominated 
form [7] and get 

A2-e5/3 

^n{^,A\^^) = ^^'{^,A';^^)-e^^A^{^,A';^^), (11) 
where we have kept only the single and double poles at t = 1 in the second term, so that 

A^(^,a2;^2) ^72^^,(a2) (^2 + ln^^ | . (12) 

In the vicinity of the pion pole one can approximate Fjj:{A^) = Ag^Mj^ / (m^ — A^). 

In estimates, shown in Fig. 2, we relied on GPDs deduced from an ansatz based on 
modeling the double distribution as a product [12] AF{y,z,A^) - 7t{y, \z\)Af{y,A^) of 
exclusive profile 7t and an inclusive parton distribution augmented to have an intrin- 
sic momentum-transfer dependence A/^(3;,A2) = rigAqj("^~"i^ (^"•^^(1 - jc)*«(l + 7^x + 
Pqy/x) with parameters fixed by the GSA forward densities [13] in A^ = limit and 
slopes = 1. 15GeV^2^ q^^^ = l.OGeV^^ chosen to fit the dipole form of the axial form 
factor with the effective mass = 0.9 GeV^. We give the cross section [at leading 
order compatible with earlier estimates [15, 6]] and transverse target-spin asymmetry 
= (2ap)/{nap). In our evaluations we set = 1 and A^s = 280 MeV for Nf = 4 
and use the tree level result for Note however that in calculations of higher- 

twist corrections via renormalons in deeply inelastic scattering in order to get the right 
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FIGURE 2. Generalized structure functions (left) in leading twist approximation (dashed) and including 
twist-four corrections (solid) as a function of xb for = -O.SGeV^ and ^ lOGeV^: (1) 5Re^ (2) 
SmJf', and (3) 10^^ • S . The photoproduction cross section in units of nbarns (middle) without (solid) 
and with (dash-dotted) power suppressed contributions for the same values of the kinematical variables. 
The transverse spin asymmetry (right) at leading order (solid) and with twist-four power effects taken into 
account for = -O.SGeV^ and = 4GeV^ (dashed) and B'^ = lOGeV^ (dash-dotted). The maximal 
value of jCB.max is set by the kinematical constaint |A^| > lAj^j^^l =M^.5t|/(l — jcb)- 



magnitude of experimental data one has to take a larger value \ Q\ ^2 — 3 [14]. The ex- 
tremely large power corrections to the absolute cross section of pion leptoproduction are 
in qualitative agreement with the earlier consideration in Ref. [16]. As we observe, how- 
ever, the renormalon model of higher-twist contributions affects in a marginal way the 
asymmetry and thus leads to the apparent conclusion of the precocious scaling in ratios 
of observables, — a fact pointed out previously in various circumstances [5, 6, 17]. 
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